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Abstract 



We derive an expression for the static dielectric constant of the colloidal 
suspensions based on the electrokinetic equations. The analysis assumes that 
the ions have the same diffusivity, and that the double layer is much thinner 
than the radius of curvature of the particles. It is shown that the dielectric 
increment of the double layer polarization mechanism is originated from the 
free energy stored in the salt concentration inhomogeniety. We also show that 
the dominant polarization charges in the theory are at the electrodes, rather 
than close to the particles. 

PACS. 77.22.Ch - Permittivity (dielectric function). 

PACS. 77.84. Nh - Liquid, emulsions, and suspensions; liquid crystals. 
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I. INTRODUCTION 



Since its inventionllrEl, the double layer polarization theory of charged colloidal suspension 
is well known for its complexity. In the theory, the dielectric function is calculated by 
analysing the out of phase (charging) electric current at finite frequencies. This approach 
makes the theory hard to grasp. Although equivalent, the out of phase current is not as 
close to the intuitive thinking as the dipole built up under the applied field. In particular 
the location of the polarization charge is not obvious from the theory. Some researchers even 
wrongly assumed that the dominant polarization charge is around the particle when trying 
to simplify the theory. 

Another consequence of using current to do calculation is that the static dielectric con- 
stant is obtained indirectly, via taking the zero frequency limit of the dynamic ac calculation. 
Normally one would like to be able to calculate the static dielectric constant statically, be- 
cause such calculation will be easier to understand than the dynamic one. However, for 
many years there is no static calculation which reproduces the known result. 

In a recent workB, Grosse and Shilov, based on the equilibrium thermodynamic argument, 
suggested that the zero frequency dielectric constant should be obtained by calculating the 
free energy stored, as the form of the non-uniform salt concentration around the particles. 
For the system of z : z electrolyte, their suggestion reads 



where Vtot is the total volume, E is the applied electric field strength. The static dielectric 
increment Ae is the static absolute dielectric constant subtracts the contribution from the 
simple electrolyte solution. The free energy stored AF, is calculated by integrating the 
free energy changes due to the counterion and colon concentration perturbations 6c^, 5c~ 
outside the double layer, where C is the equilibrium ion concentration in the bulk, k, T are 
the Boltzmann constant and the temperature respectively. In fact, the ions concentrations 
are the same outside the double layer, therefore the free energy is stored in the form of 
non-uniform salt concentration distribution. 
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Grosse and Shilov have checked that eq.(|T|) is obeyed for the well known solution of the 
dilute charged sphere suspension. They further used eq.([l|) to calculate the static dielectric 
increment of dilute ellipsoids, which is not easy if the original calculation method is used. 

In this work, we prove eq.(|l|) (with some corrections from the double layer region) based 
on the electrokinetic equations, so that its general validity is established on the same basis 
as the double layer polarization theory. We consider only the simple case where the two ions 
have the same mobility, but put no restriction on the C potential value. By analysing the 
charging current, we also show explicitly that the polarization charges are close to the two 
electrodes. Below, we summarize the thin double layer theory in section II. The charging 
electric current density and the dipole are analysed in section III. Eq.(|I]) is derived in section 
IV. A short summary concludes the paper. 



II. THE THIN DOUBLE LAYER THEORY 



The standard electrokinetic theorycTu'a describes the ion dynamics (5c^) by the Smolu- 
chowski equations, the electrostatic potential (0) follows the Poisson's equation, and the 
fluid velocity obeys the Navier-Stokes equation with an electric body force given by the 
product of the local charge density and the electric field. 

The analysis of the dynamic equation proceeds differently within and outside the double 
layer (the bulk), where different approximations been appropriateBi'i. The bulk equations 
to be used later are summarized in subsection [11 A] below (which apply for the z : z elec- 



trolyte with the equal ion mobility). The partial solutions from the two regions are matched 
smoothly at the boundary of the double layer to give the full solution. 

In the case of a thin double layer, where the radius of curvature of the particle is much 
larger than the double layer thickness, the double layer region can be approximated as a 
flat double layer. The solution within the fiat double layer is known and the matching 
conditions act effectively as the boundary conditions for the quantities outside the double 
layer. In subsection [II B| below we summaries the boundary conditions to be used. 



3 



A. The Bulk Equations 



It is well known that, in the case of z : z electrolyte with the equal ion mobility, the 
solution becomes charge neutral a few Debye lengths away from any interface0. The salt 
perturbation does not carry charge, therefore uncouples with the potential. The unperturbed 
salt has a uniform distribution so the flow does not couple with the salt perturbation by 
convection (to the first order perturbation). The bulk solution is charge neutral to begin 
with, so the potential does not produce force to couple with the flow. In the bulk one only 
needs to solve the Laplace equation 

= (2) 

for the potential, and the simple diffusion equation 

dt6c = DV'^dc (3) 

for the salt perturbation 6c = {6c^ + 6c^)/2. The diffusivity D is the same for both ions. 
Since the flow velocity couples passively to other quantities (via the boundary condition at 
the double layer), we do not need to calculate it for analysing the dielectric response. 
The conduction current density i in the bulk is local and Ohmic, obeying 

i = KE (4) 

The conductivity is given hy K = 2Dq'^C/kT, where q is the counterion charge. Note that 
eq.(^) is justified for the electrolyte with the equal ion mobilities, since the perturbed salt 
is charge neutral. For systems with different ion mobilities, there is also a (small) current 
density controlled by the salt, which is not proportional to the local electric field. The fact 
that ions in general have different mobilities complicate the analysis very much (all three 
eqs.(0)(^(§) will be modified). The correction in the dielectric response!, however, seems 
perturbative. We will therefore focus on the simplest system in this work. 



4 



B. The Double Layer Boundary Conditions 

Define the perturbed surface cfiarge density as cr = g ^q^'^{5c^ — 6c^)drn wfiere r„ is tlie 
coordinate in tlie direction normal to tlie interface (which is r„ = 0). By analysing the ion 
conservations within the thin double layer, and assuming that no ion fluxes penetrate into 
the particles, one obtains the conservation equation for the surface charge density 

dta = KVn(p + KkiiV\(l) + Kki2VlC (5) 

where C = kT5c/qC. The Debye length is = {2q^C / ekTy/'^ where e is the dielectric 
constant of the solvent. Note that a is not sensitive to the exact value of the cut off a, 
whose plausible value is between 1 to 10. The surface conductances are expressed as the 
bulk conductance K times the relative conductances fcjjS, to make some formulae neater. 

Using a similar analysis, the perturbed surface excess ion density 7 = /q"^'" [5c^ + 5c~ — 
25c(r„ = a/K)]drn obeys 

qda = KVnC + Kk2iVl4> + Kk22VlC (6) 

Here 7 is also not sensitive to the exact value of the cut off. 

The surface conductances kij have been calculatedBi, which include the effect of the 
convection current inside the double layer. There is also an Onsager relation ku = ^21 as 
can be checked in the microscopic expressions!. There are also other boundary conditions for 
the flow field, and constitutive equations for a, and 7. However we will only need eqs.(||)(|^) 
in this work. 

Note that in the double layer polarization theory, the slow polarization of the double 
layer is controlled by the salt diffusion. Therefore the time derivative terms dtcr, dt'~f in 
eqs.(^)(^ are often small, hence neglected. In this paper the time derivatives are kept, so 
that the interface polarization0 (Maxwell- Wagner theory) is included. 
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III. CHARGING CURRENT DENSITY AND THE DIPOLE 



In this section, we decompose the total current density into the charging current density, 
which stores free energy, and the conducting current density, which is dissipative. We 
also show explicitly that the large dipole comes from the polarization charge close to the 
electrodes, rather than from the charge close to the particles. 



A. Averaged current density 

Before we decompose the total current density, we first express the averaged current 
density as did by De Lacey and White!. The final expression has the advantage that it is 
valid independent of the specific current density constitutive relation close to or within the 
charge particles. 



VtoSa) = f iadV+ [ iadV 

= f iadV + f \/biraib)dV - f raV-idV 

Jv,„, JVp JVp 



f KEJV) - 


f /" rJ-dA+ f 


IVsol J 





= (k [ EJV-K I Ejv] + I rJ-dA+ [ r^dtPfdV 

\ JVtot JVp J JSp JVp 

= K f EJV + K f {<f)6at - r„V50) dA^ + f r^dtPfdV (7) 

JVtot J Sp JVp 

We define a two- components surface Se, which is inside the solution and very close to the 
two electrodes (see Fig. 1). The surface encloses the total space Vtot, which is further 
partitioned by the surface Sp into Vtot = Vp + Vsoi- Sp is a (multi-components) mathematical 
closed surface which includes the charged particles and their double layer (Fig. 1). The space 
Vp is enclosed by the surface Sp. The volume Vgoi is the space between the two surfaces Sp 
and Se- The differential area vector dA is pointing into Vsoi for both surfaces Se and Sp. 

In the above derivations, the space integral is decomposed into two term. The integral 
on Vp is transformed into two term by the identity ia = Vb{raib) — faVbh- To get the third 
line we have assumed that the current in the volume Vsoi obeys the Ohm's law i = i^E. 



The second term is obtained by the Gauss' theorem. In deriving the last term, we have 
not assumed any constitutive relation for the current density in Vp. The only property used 
inside Vp is the charge conservation dtPf = —V ■ i where p/ is the free charge density. In 
the fourth line the integral on Vsoi is written as Vtot ~~ Vp. In the last line, the relation 
Ea = — Va0 and the Gauss' theorem are used. We also assumed that the Ohm's law applies 
on the surface Sp. 

B. Current density decomposition 

For each instance, we define the conducting current density i° as the current density of a 
reference system, which is at a steady state where the averaged electric field is the same as 
the original system at that time. The charging current density 6i is the difference between 
the total current density and the conducting current density. Given a time evolution of 
the applied field, for example an adiabatic charging process, the conduction current density 
as defined does not change sign under the time reversal operation. Whereas the charging 
current density changes sign under the time reversal. 

The reference system has the fields 0°, C° which both obey the Laplace equations, and 
the steady state boundary conditions 

irVn0° + KkuV\(tP + Kki2V\C^ = (8) 
KV„C° + irA;2iVi0° + Kk22VlC^ = (9) 

The analysis similar to (0) on the conducting current i° gives 

Vtotii^) = K I YPdV + K I f0°l - rV0°) • dA (10) 

JVtot JSp ^ ' 

where 1 is the identity matrix. Equation (|TD|) has one term less than eq.(|^) because V-i° = 0. 

Define the charging current density as 5i = i — i*^. Noticing that /y^^^ Ec/V^ = /y^^^ Y^^dV 
(by the definition of 0°), we subtract eq.(0) from eq.(^ to cancel their first terms to give 

Vtot{^\) = K f (501 - rV50) ■dA+ f rdtPfdV (11) 



where 50 = — 0° and 6C = C — obey the Laplace and diffusion equations respectively. 
Their boundary conditions read 

dta = KVJ(j) + KknVl6(f) + KkuVldC (12) 
qda = KVJC + Kk2iV\5(t) + Kk22V\5C (13) 

It is the first term K J^^ {S(j)l — rV(50) ■ dA in eq.(0) which gives the main contribution 
in the salt polarization theoryi. 

C. Charging current density and the electric dipole 

It is instructive to examine the dipole produced by the charging current density {Si). 
Vtot{Sia) = / SiadV 

JVtot 



f Vb{ra6tt)dV - f r^V-idV 

J Vtot J Vtot 

- / TaSi ■dA+ TadtpfdV 

J Se JVtot 



where we have used the identity ia = Vb(ra4) — raVbh to get the second line. In the last 
step, Gauss's theorem and the charge conservation equation are used. In the last term, we 
can restrict the range of integration to Vp because the charge density pf vanishes within Vsoi 
for the symmetric electrolyte. The result reads 

dtPtot = Vtot{Si) = - f r5i-dA+ [ rdtPfdV (14) 

where 'Ptot is the dipole of the system. The first term — /^.^ rS'i ■ dA is the dipole increment 
due to the free charges flow into the electrode surface. The term /^^ rdtpdV is the time 
derivative of the dipole close to or inside the particles. Compare eq.(0) with eq.(|l4|), we 
conclude that the large polarization in the double layer polarization theory consists of the 
dipole built up by the (free) ions accumulated at the boundary of the electrodes, rather than 
the dipole around or within the particles. 
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IV. STATIC DIELECTRIC CONSTANT AND THE FREE ENERGY STORED 



We now establish our main goal, the reversible work done at the charging process is due 
to the free energy stored as the form of the non-uniform salt concentration. Our result also 
contains the additional contributions within the double layer regions. The charging power 
iJ 

V = f (j)6i-dA (15) 

'J Se 

Note that Se has two components close to the two electrodes. Equation ([ISD is identical to 
the charging current times the voltage difference between the two electrodes. 

Because (E) = (E°) and both 0, 0° are constant on each electrodes, the reference 
potential 0° is the same as the true potential at the electrode surface Sg. Therefore can be 
replaced by 0° in eq.(|T5|). In the bulk solution, the current is simply given as 6i = —KVdip. 
Using that V^0° = V^(50 = within Vsou one integration by parts to get the identity 

/ 0°V50 ■dA+ [ (tPV5(t) ■dA = - [ V ■ f0°V50) dV = - [ V0° ■ V6(j)dV 

JSe JSp JVsol ^ ' JVsol 

= - [ V ■ f 50V0°) dV = [ 6(f)V(f)° ■dA+ [ 5(t)V(tP ■ dA 

The last term is actually zero because 50 = at the electrodes 5*6. Use the above identity, 
we obtain 

V = -K I (lPV5(j) ■dA = K I {4PV5(I) - 5(j)V(f)^) ■ dA (16) 

We first rearrange the term / 0°V50 ■ dA so that the driving force of the surface current 
0° + (^12/^11)^° appears (see the last two terms of eq.(H)). 

/ 0OV50-rfA=/ f0O + :^C°W50-rfA-^ / C°V50-rfA 

JSp JSp \ kii J kii JSp 

In the first term, V(50 ■ dA = V„50(iA, and we eliminate V„50 by eq.(0), to get 

4/ (<P' + ^Ad,adA-h, [ (<P' + ^Avl(5<P+^5c]dA-^ [ C°V50-rfA 
K JSp \ kii J JSp \ kii J \ kii ) fell JSp 



JSp \ kii ) Jsp \ ku J \ fell J fell Jsp 



V64> ■ dA 



where (surface) integration by parts is used twice in the second term. (Note that Sp is a 
close surface, and has no boundary.) Finally we use eq.(|D to eliminate 0° + (A;i2/fcii)C° in 
the second term to get 



The first term can be further transformed by using eqs.(D(0) to eliminate V„0°, dtcr — 
KVn^cj) respectively, as 



where Jg^ SCV^C^dA = Jg^C^V^SCdA is used to cancel terms. In the last step, (PD (p!3D are 
used to eliminate V^^", V\6(f) respectively. 

We would like to enlarge the integration area of the first two terms to include Se, so 
that the surface integral on Sg and Sp can be transformed to the bulk integral on Vsoi by 
the Green's identity. The salt concentration C° contributed from each particle decays as 
or faster, because the monopole strength vanishes § dA ■ VC° = (we use eq.(P), and 
integration by parts). The quantity 6C decays no slower than C°. Therefore the integral on 
the electrode surface Se scales less than the system size. Whereas the original integrals on 
Sp are extensive, and hence dominant. For a large enough system (thermodynamic limit), 
we freely enlarge the integration area of the first two terms as Sp + Se with negligible error, 
to derive 



1 / f / + d,adA + f f 50 + ^6C] Vn<P'dA - ^ / C'VJHA 

K JSp \ kii j JSp \ kii J kii JSp 



The above manipulations lead to 





k 






rdt6CdV 
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where we have used the Green's identity, V^C° = 0, and dt5C = DV^5C. 

Finally, put the above derivations together, and use the Onsager relation ki2 = k2i, we 
obtain the result 

K 



— I C^dtdCdV + I Cydt-idA + / (l)^dtadA 
f ^dt6c+ I ^dt^dA+ I <pdtadA (17) 



where in the last step, we have assumed that the charging rate is very small, to neglect the 
terms quadratic in the charging rate. 

In the linear response regime, the perturbations 7, a, 6c, and are all proportional. We 
can integrate eq.(|T^ to obtain 



Vtot— 



-E' = AF = / — ^rfl^ + / ^idA + / ^adA 
JVsol JSp 2C JSp 2 



Equation dlSf ) shows that the external electric work done of the charging process is actually 
saved in the forms of the non-uniform salt concentration in the bulk (first term) and inside 
the double layer (second term), as well as the electric energy stored in the polarized double 
layer (last term). The first term is as proposed by Grosse and Shilovi. The second term is 
of the similar physical origin as the first term but its magnitude is much smaller for large 
particles with thin double layers. The dynamic calculationsS'i do not have this term because 
they neglect the (small) term dtj in eq.(^ so that their static result is less complete as in 
eq. (|T^) . The last term accounts for the energy stored in the interface polarization^. It is 
not included in the dynamic calculation&i because the term dta in eq.(|^) has been neglected. 
The exact values of a and 7 depend on the quantities like the dielectric constant and the 
conductivity of the particles, but the general forms of the last two terms in eq. (|TBp remain 
the same. 
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V. SUMMARY 



Using the electrokinetic equations, we have shown in eg. ([T8|) that the work done on the 
charging current is stored in (i) the free energy of the uneven salt concentration in the bulk 
solution, (ii) the free energy of the double layer excess ion density, and (iii) the polarized 
double layer charge. 

We also show in eqs.(|Tl|)(|l^ that, the dipole of the free charge close to the electrode 
surface, gives the dominant dielectric increment of the double layer polarization theory. 
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FIGURE CAPTION 



Figure 1. Two surfaces Sp and Sg are shown. The space Vsoi is bounded between the 
two surfaces. The space Vp is within the surface Sp, which includes the charged particles 
and their double layers. 
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